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ABSTRACT 



The Far Ultraviolet Spectroscopic Explorer (FUSE) has surveyed a large sam- 
ple (> 100) of active galactic nuclei in the low-redshift universe (z < 1). Its 
response at short wavelengths makes it possible to measure directly the far ultra- 
violet spectral properties of quasistellar objects (QSOs) and Seyfert 1 galaxies at 
z < 0.3. Using archival FUSE spectra, we form a composite extreme ultraviolet 
(EUV) spectrum of QSOs at z < 0.67. After consideration of many possible 
sources of systematic error in our analysis, we find that the spectral slope of the 
FUSE composite spectrum, a = — 0.56tg!| f° r ^ v ^ v<X i ^ s significantly harder 
than the EUV (A < 1200 A) portion of the composite spectrum of QSOs with 
z > 0.33 formed from archival Hubble Space Telescope spectra, a = —1.76 ±0.12. 
We identify several prominent emission lines in the FUSE composite and find 
that the high-ionization O VI and Ne vm emission lines are enhanced relative 
to the HST composite. Power law continuum fits to the individual FUSE AGN 
spectra reveal a correlation between EUV spectral slope and AGN luminosity in 
the FUSE and FUSE+HST samples in the sense that lower luminosity AGNs 
show harder spectral slopes. We find an anticorrelation between the hardness 
of the EUV spectral slope and AGN black hole mass, using estimates of this 
quantity found in the literature. We interpret these results in the context of the 
well-known anticorrelation between AGN luminosity and emission line strength, 
the Baldwin effect, given that the median luminosity of the FUSE AGN sample 
is an order of magnitude lower than that of the HST sample. 
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1. Introduction 

The ubiquity with which QSOs display 
spectral properties such as power law con- 
tinua and broad emission lines over wide 
ranges in luminosity and redshift has led to 
the use of composite spectra to study their 
global properties. In this paper, we present 
a composite far, or extreme, ultraviolet 
(EUV) spectrum of low-redshift AGNs. In- 
formation about the continuum in the rest- 
frame ultraviolet is particularly critical for 
understanding the formation of the emis- 
sion lines, for characterizing the UV bump 
in QSO spectral energy distributions, and 
for determining the ionization state of the 
intergalactic medium (IGM). Composite 
QSO spectra covering the rest-frame ul- 
traviolet have been constructed for AGNs 
with 0.33 < z < 3.6 from HST (Zheng et 
al. 1997; Telfer et al. 2002, Z97 and T02 
hereafter) , and at z > 2 from ground-based 
samples like the Large Bright Quasar Sur- 
vey (Francis et al. 1991), the Sloan Digi- 
tal Sky Survey (SDSS, Vanden Berk et al. 
2001), and the First Bright Quasar Survey 
(Brotherton et al. 2001). Z97 and T02 re- 
ported that the HST composite spectrum, 
which covers rest wavelengths ~350-3000 
A, shows a spectral break in the power 
law continuum at 1050-1300 A in the sense 
that the EUV spectral shape blueward of 
the break is softer than the slope red- 
ward of the break, in the near ultravio- 
let (NUV). In this composite, they iden- 
tify several prominent far-UV emission fea- 
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tures including lines due to Lya blended 
with the N v doublet, Ljj3 blended with 
the O VI doublet, and Ne VIII. 

The bandpass of the Far Ultraviolet 
Spectroscopic Explorer (FUSE) (Moos et 
al. 2000; Sahnow et al. 2000), 905-1187 
A, allows us to examine the EUV prop- 
erties of local AGNs. Therefore, we can 
study the same rest-frame EUV wave- 
length region covered by the HST com- 
posite spectra, for AGNs with redshifts 
less than 0.33. The low redshifts of these 
FUSE AGNs ensure that although the 
FUSE aperture limits it to observing rel- 
atively bright AGNs, our sample contains 
a substantial fraction of intrinsically low- 
luminosity (log(AL A [erg s -1 ]) at 1100 A 
= log(ALnoo) < 45) AGNs. The FUSE 
sample combined with the HST sample 
yields a sample of AGNs with a spread of 
nearly five orders of magnitude in luminos- 
ity from which we can investigate trends 
in EUV spectral shape with luminosity in- 
voked to account for the anticorrelation 
between AGN luminosity and emission line 
equivalent width, the Baldwin effect (Bald- 
win 1977; Netzer, Laor, & Gondhalekar 
1992; Green 1996, 1998; Wang, Lou, & 
Zhou 1998; Dietrich et al. 2002), and the 
dependence of the strength of the Bald- 
win effect on the ionization potential of the 
emitting ion (Zheng, Fang, & Binette 1992; 
Zheng & Malkan 1993; Zheng, Kriss, & 
Davidsen 1995; Espey & Andreadis 1999; 
Dietrich et al. 2002; Kuraszkiewicz et al. 
2002; Shang et al. 2003). 

The AGNs in our FUSE archival sam- 
ple all have redshifts less than 0.7, afford- 
ing the advantage that the determination 
of the mean EUV spectral index requires a 
less significant correction for IGM absorp- 
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tion than was required for the HST sam- 
ple. After consideration of various possible 
systematic effects on the analysis, we com- 
pare the spectral shape of the composite 
and the strength of the emission lines to 
those found in the EUV spectrum of AGNs 
with z > 0.33 compiled from HST data 
by T02. We also fit power law continua 
to each individual AGN spectrum in the 
FUSE sample and examine the results for 
correlations of the spectral slope with red- 
shift and luminosity. Finally, we compile 
estimates of black hole mass for several of 
the sample AGNs to test for an anticor- 
relation between this quantity and spec- 
tral slope. Such a correlation is a predic- 
tion of an evolutionary model in which the 
central black holes of AGNs accrete mass 
over time and the peak of the UV bump in 
the spectral energy distribution shifts to 
longer wavelengths, resulting in softer ion- 
izing continua in the UV and soft X-rays 
for higher mass/luminosity AGNs (Wan- 
del 1999a,b). In such a scenario we would 
also expect to find a correlation between 
UV spectral slope and the accretion disk 
temperature, as estimated from the AGN 
luminosity and black hole mass, and we in- 
vestigate this with the FUSE data as well. 

2. The Sample 

From the FUSE archives, we down- 
loaded the 165 spectra of AGNs with z < 1 
that were public as of 2002 November. 
We processed the raw data using stan- 
dard FUSE calibration pipelines (see Sah- 
now et al. 2000) to extract the spectra, to 
perform background subtraction using up- 
dated background models and subtraction 
algorithms, and to perform wavelength and 
flux calibrations. We also implement a 



correction for the "worm", a dark stripe 
running in the dispersion direction on one 
of the FUSE detector segments. See the 
FUSE web pages" 5 for details about this 
feature. 

Following a procedure similar to that 
of T02, we excluded spectra of broad- 
absorption-line quasars and spectra with 
signal-to-noise ratios (S/N) < 1 over large 
portions. We also exclude spectra of AGNs 
with strong narrow emission lines, e.g. 
NGC 4151, strong stellar features in their 
continua, e.g. NGC 7496, or strong inter- 
stellar molecular hydrogen absorption, e.g. 
1H 2107-097. A total of 128 spectra of 85 
AGNs, all with z < 0.67, meet the criteria 
for inclusion in the sample. The 85 AGNs 
in the FUSE sample are listed in Table 1, 
along with their redshifts and Galactic red- 
dening values from Schlegel, Finkbeiner, & 
Davis (1998). 

3. Composite Spectrum Construc- 
tion 

We follow the same procedure as T02 
for the construction of the composite spec- 
trum. The reader should consult that pa- 
per for a detailed discussion of the cor- 
rections applied to each sample spectrum, 
which we summarize here. 

1. We correct the spectrum for Galac- 
tic extinction using the Cardelli, 
Clayton, & Mathis (1989) extinc- 
tion curve, E(B — V) values listed 
in Table 1 (Schlegel et al. 1998), and 
R v = 3.1. 

2. We exclude wavelength regions af- 
fected by interstellar absorption lines. 

5 http: / / fuse.pha.jhu.edu 
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3. We correct for Lyman limit absorp- 
tion using r LLS = ln(F + /F_), if 
the S/N below the Lyman break is 
greater than one, where F + is the 
median flux in selected windows red- 
ward of the Lyman break and F_ 
is the median flux blueward of the 
break. 

4. We apply a statistical correction for 
the line of sight absorption due to 
the Lya forest and the Lyman val- 
ley (M0ller & Jakobsen 1990). 

5. We shift the AGN spectrum to the 
rest frame. 

6. We resample the spectrum to com- 
mon 1 A bins. 

The lower redshifts of our sample AGNs 
compared with the HST sample of T02 
compels us to use different parameters to 
correct for Lja forest absorption men- 
tioned above, so we discuss them briefly 
here. Like T02, we use the distribution of 
absorbers given by 

^«d + W (i) 

We account for column densities in the 
range 12.2 < logiV < 16.7. For the col- 
umn density distribution parameter, we 
use the result found by Dave & Tripp 
(2001) from iJST/Space Telescope Imag- 
ing Spectrograph (STIS) echelle spectra of 
two QSOs at z ~ 0.3, (3 = 2.0 for 12.2 < 
logiV < 14.4. For 14.4 < logiV < 16.7, we 
use (3 = 1.35 from the //iST/Goddard High 
Resolution Spectrograph study by Penton 
et al. (2000). For the redshift distribution 
parameter, we use 7 = 0.15 (Weymann et 
al. 1998). We normalize the distribution 



of absorbers in Equation 1 by 1.34 x 10~ n 
cm 2 at logiV = 13 and z = 0.17 and as- 
sume a Doppler parameter of 21 km s _1 
(Dave & Tripp 2001). 

We combine the sample spectra using 
the bootstrap technique described by T02. 
To summarize briefly, we begin the boot- 
strap procedure at the central portion of 
the output composite, specifically the re- 
gion between 850 and 950 A. We then 
include spectra that fall at longer wave- 
lengths in sorted order to longer wave- 
lengths. Finally, we include those at 
shorter wavelengths in sorted order to 
shorter wavelengths. The overall compos- 
ite is renormalized at each step. Figures 1 
and 2 show the number of spectra con- 
tributing to the final composite in each 
wavelength bin, and the S/N per wave- 
length bin in the final composite, respec- 
tively. The non-smooth appearance of 
these two histograms is caused by the large 
number of discrete wavelength regions sur- 
rounding interstellar absorption lines we 
omitted from the calculation. 

The final FUSE composite spectrum, 
shown in the top panel of Figure 3, cov- 
ers the rest wavelength rang e 630-1155 A. 
We used the IRAF task specfit (Kriss 1994) 
to fit a power law of the form F v oc v a 
to the continuum of the composite using 
wavelength regions free of emission lines: 
630-750, 800-820, 850-900, 1095-1100, and 
1135-1150 A (T02). The best-fit power 
law index is a = —0.56, and we show 
this fit with the dashed line in Figure 3. 
We compare this index to the values of 
«euv from continuum fits to the HST com- 
posite for wavelengths > 500 A in T02. 
In the top panel of Figure 3, we also 
show the HST composite from T02, and 
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in the bottom panel of the figure, we show 
the ratio of the FUSE and HST compos- 
ites. The FUSE composite is significantly 
harder and shows enhancement of several 
emission lines, which will be discussed fur- 
ther in Section 6. 

4. Uncertainties 

Because the S/N in the FUSE compos- 
ite is high (Figure 2), the statistical error 
in the continuum fit is small. The largest 
sources of error are likely to come from 
cosmic variance in the spectral shapes of 
individual AGNs and from systematic er- 
rors in the method we have used to correct 
the sample spectra for Galactic and inter- 
galactic absorption. We estimate the error 
arising from the range of spectral shapes 
of the individual AGNs that constitute our 
FUSE sample by creating 1000 bootstrap 
samples with replacement from the origi- 
nal sample. From the bootstrap samples, 
we find a standard deviation in a of 0.11. 

We now explore how a number of pos- 
sible systematic errors within our analy- 
sis would affect the results. The boot- 
strap combination technique is robust, so 
varying it has little effect on the resulting 
composite. If we perform the bootstrap 
from long wavelengths to short, a changes 
by +0.03. Normalizing the spectrum at 
different wavelength intervals, 1050-1150 
A (T02) or 800-900 A versus the fiducial 
range 850-950 A changes a by an amount 
less than the statistical error in the power 
law fit from specfit, +0.01/ - 0.003. 

We find that the results are sensitive to 
the extinction correction in terms of both 
the individual values of E(B — V) and 
the adopted ratio of total-to-selective ex- 
tinction, R v . Changing all individual val- 



ues of E{B — V) by ±lcr, where we esti- 
mate la = 0.1QE(B - V) (Schlegel et al. 
1998) changes a by ±0.16. For our fidu- 
cial correction for Galactic reddening, we 
used the average value of R v for sightlines 
through the diffuse interstellar medium of 
the Milky Way, 3.1 (Cardelli, Clayton, & 
Mathis 1989). However, for individual 
lines of sight, it may vary as widely as 
R v = 2.5 to R v = 5.5 (Clayton & Cardelli 
1988). If instead of Rv = 3.1, we use in- 
stead Rv = 2.8 or Ry = 4.0 in the extinc- 
tion law, the spectral index of the compos- 
ite spectrum changes by +0.06 and —0.19, 
respectively. 

The composite is also sensitive to the 
value of the column density distribution 
parameter, f3. Because the result of Dave 
& Tripp (2001) is larger than typically 
inferred for f3 from lower resolution data 
(Penton et al. 2000; Dobrzycki et al. 2002) 
or at high redshift (Kim, Cristiani, & 
D'Odorico 2001), we reduce it from the 
fiducial value of 2.0 to 1.5. This increases 
a by 0.3. The composite is fairly insensi- 
tive to the redshift distribution parameter, 
7, which we varied from 0.15 (Weymann et 
al. 1998) to 0.54 (Dobrzycki et al. 2002). 
Raising and lowering the upper column 
density limits for the Lya forest correction 
to logiV = 16.9 and logiV = 16.5 changes 
a by +0.06 and —0.04, respectively. 

The relatively strong effect that varying 
j3 has on the results is somewhat surpris- 
ing, since at the low redshifts considered 
here, the number of Lya forest absorbers in 
any given spectrum is small. At the far ul- 
traviolet wavelengths relevant to this com- 
posite spectrum, we have corrected also for 
line blanketing in the Lyman series and for 
Lyman continuum absorption, the Lyman 
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valley (M0ller & Jakobsen 1990), using the 
same Lja forest parameters as discussed 
above. The effect that changing these pa- 
rameters has on the Lyman valley correc- 
tion for a z = 0.1 AGN is shown in Fig- 
ure 4. We also show here the correction 
calculated using the Lyct forest parameters 
of T02. This correction is less than 1% over 
the whole spectrum because, relative to the 
shallow power law number distribution ob- 
served at low redshift, T02 used 7 = 2.46 
appropriate for the high redshifts of their 
HST AGN sample. Extrapolating that 
power law to z — 0.1 with the same nor- 
malization results in a factor of ~50 under- 
estimate in the number of Lya absorbers 
at z — 0.1. In this figure, we can see the 
effect that changing 7 and (3 has on the 
correction used for the FUSE data. Chang- 
ing 7 from 0.15 to 0.54 has relatively little 
effect, while increasing (3 from 2.0 to 1.5 
increases the correction substantially, to 
>20%. It should be noted, however, that 
in any given AGN spectrum, no identifiable 
Lyman valley trough is observed. Deriving 
a composite spectrum with no correction 
for the Lja forest and Lyman valley re- 
sults in a = —0.75, a marginally softer in- 
dex than quoted above, though still signifi- 
cantly harder than the EUV spectral index 
derived from the HST data by T02. 

Adding all the systematics discussed in 
this section in quadrature, we estimate the 
total uncertainty in a to be (+0.38, —0.28). 
The spectral index from the HST compos- 
ite, a = —1.76 ± 0.12, differs from the 
FUSE value, a = —0.56, by ~3 times this 
amount. 



5. Redshift, Luminosity, and Red- 
dening Subsamples 

Here we examine trends in the FUSE 
AGN sample with redshift, luminosity, 
and Galactic reddening using composites 
formed from various subsamples of the full 
AGN sample. The subsamples discussed 
in this section are summarized in Table 2. 

First, we calculate the luminosity of 
each FUSE AGN using the prescription de- 
scribed by T02: 

L A (1100 A) = F c (1100 A)^i^i47rDl, 

(2) 

where F c is the flux in the composite, F s 
is the flux in the individual spectrum, and 
the sum is performed over the spectral re- 
gions where the individual AGN spectrum 
overlaps with the composite. The term 
Dl is the luminosity distance to the AGN. 
These luminosities are listed in Table 1. 
We assume f2 = 1 and H = 60 km s _1 
Mpc" 1 , as T02 did. We show the distri- 
bution in redshift and luminosity of the 
AGNs in the FUSE sample and in the HST 
sample of T02 in Figure 5 with the me- 
dian redshift and luminosity of the FUSE 
sample, z = 0.10 and log(AL 110 o) = 45.0, 
marked by the vertical and horizontal lines. 
Figure 6 shows these distributions in his- 
tograms. The redshift histogram illus- 
trates that the FUSE sample is highly con- 
centrated around the median redshift. The 
luminosity histogram shows the large span 
of UV continuum luminosity in the com- 
bined FUSE + HST sample. 

Dividing the sample in redshift at z — 
0.10 and constructing composite spectra 
from the two redshift subsamples, we find 
the same value of a for both low-redshift 
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and high-redshift composites, a = —0.74. 
With this redshift cut, the total compos- 
ite spectrum is dominated by low-redshift 
AGNs at A > 950 A and by high redshift 
AGNs at A < 950 A. Figure 7 shows the 
composites made from all AGNs, and from 
high-z and low-z AGNs. 

Dividing the sample in luminosity at 
log(ALiioo) = 45.0 gives a = 0.19 for 
low-luminosity AGNs and a = —0.84 for 
high-luminosity AGNs. Figure 8 shows the 
composites made from all AGNs, and from 
high-luminosity and low-luminosity AGNs. 
Similar to the low- and high- redshift sub- 
samples, the high- and low- luminosity 
subsamples cover different spectral ranges. 
Only a small spectral region is common to 
both subsamples: 866-901 A. The differ- 
ence between the continuum slopes of the 
low-luminosity and high-luminosity com- 
posites may come from the limited limited 
spectral range used in the fit in each case. 
However, the difference between these two 
slopes is 2.7o" + , where a + is the overall un- 
certainty in the spectral index of the FUSE 
composite discussed in Section 4. The low- 
luminosity composite has a harder spectral 
shape, and Figure 8 illustrates that it also 
shows enhanced O Vl/Ly/3 emission. 

The AGNs in the FUSE sample lie along 
sightlines with Galactic reddening values 
in the range 0.006 < E(B - V) < 0.182. 
The median value is E(B-V) = 0.032. We 
show the distribution of reddening values 
in Figure 9. Because we found the over- 
all FUSE composite to be sensitive to the 
reddening correction, as discussed in Sec- 
tion 4, we consider whether a bias is in- 
troduced into the composite spectrum by 
over- or under-correcting for Galactic red- 
dening by dividing the sample into two 



subsamples, AGNs with E(B-V) > 0.032 
and those with E(B — V) < 0.032. For 
these two subsamples, we calculate a com- 
posite spectrum and find a = —0.21 for the 
high-E(B — V) composite, and a = —0.50 
for the \aw-E(B — V) composite. This dif- 
ference, in the sense that the high- E (B — 
V) composite has a harder spectrum, may 
indicate that we are over-correcting for 
Galactic reddening. However, the differ- 
ence is only marginally significant, ~ ler~ 
where a~ is the overall uncertainty for the 
composite based on the total FUSE sam- 
ple quoted above in Section 4. The distri- 
bution in redshift and luminosity of these 
two subsamples is shown in Figure 10. For 
the high-E(B — V) subsample, the median 
redshift and luminosity are z = 0.06 and 
log(AL 110 o) = 44.5, respectively, and the 
median redshift and luminosity are z = 
0.16 and log(ALnoo) = 45.5 for the low- 
E(B — V) subsample. Thus, the differ- 
ence in spectral index between the high- 
and low-E^-B — V) subsamples may also 
reflect an underlying trend toward harder 
EUV spectra in lower-luminosity AGNs. 

6. Emission Lines 

We fit emission lines to the FUSE com- 
posite spectrum using specfit. We show 
these lines in Figure 11, and their param- 
eters are listed in Table 3. The fluxes are 
normalized such that the O VI emission line 
has a flux of 10. We identify Ne Vlll at 774 
A, O ill at 831 A, a broad feature at 944 A, 
attributed to the H I Lyman series with a 
possible contribution from S VI, Ly7 at 973 
A, C ill and N ill at 976 and 991 A, respec- 
tively, Ly/3 plus O VI at 1026 and 1033 A, a 
S IV doublet at 1062 and 1073 A, and a fea- 
ture at 1084 A, possibly due to Nil, Hell, 
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and/or An (Z97; T02). A dip is present in 
the FUSE composite spectrum blueward of 
the Ne VIII line at ~730 A. This feature is 
visible also in the HST composite from T02 
in Figure 3. 

The bottom panel of Figure 3 shows 
that LjP + O VI and Ne VIII are strongly 
enhanced relative to the HST sample. The 
equivalent widths of these lines are ~70% 
and ~90% larger in the FUSE composite. 
These spectral enhancements are the best 
piece of evidence that the bluer continuum 
of the FUSE composite is not due to a sys- 
tematic error in the corrections for Galac- 
tic reddening and/or IGM absorption but 
rather due to real physical differences in 
AGN properties. We shall discuss this fur- 
ther in Section 8. 

7. Spectral Fits to Individual AGNs 

We fit a power law to each AGN EUV 
continuum individually for the 85 FUSE 
AGN spectra. The spectral slopes, a\, and 
extinction-corrected normalizations, F\ at 
1000 A, for the fits to F x oc X~ ax are listed 
in columns 5 and 6 of Table 1. In the fol- 
lowing discussion, we refer to the spectral 
slopes of the individual AGN as defined 
above, F v oc u a , i.e. a = oc\ — 2. Figure 12 
shows the distribution of the EUV spectral 
slopes in the FUSE sample and in the com- 
bined sample from FUSE and HST, with 
median values for each sample marked. We 
will discuss this combined sample further 
in Section 7.3.1. 

7.1. Correlations of Spectral Index 
with Reddening and Flux 

We return once more, briefly, to the is- 
sue of bias introduced by our Galactic red- 



dening correction by examining the indi- 
vidual AGN spectral slopes for any trend 
with E(B - V). A plot of a versus E(B - 
V) is shown in Figure 13. From a least- 
squares linear fit to these points, we find 
a slope of 17.4 ± 0.5, but the Spearman 
rank-order test indicates that the correla- 
tion is not significant (70%). This result 
gives us further confidence that our over- 
all composite is not significantly biased by 
errors in our correction of the individual 
FUSE spectra for Galactic reddening. 

We now investigate whether our results 
may be influenced by systematic errors in 
the background subtraction we performed 
on the FUSE data. This is a particu- 
lar concern for faint AGNs, especially in 
the short wavelength regions of the spec- 
tra, covered by the SiC channels of the 
FUSE spectrograph, which have lower ef- 
fective area than the LiF channels by a 
factor of ~3. If we have systematically 
over- or under- subtracted scattered light 
in the SiC channels, we expect a correla- 
tion between a and the observed flux. We 
show the individual spectral slopes versus 
the flux at 1000 A from the power law 
continuum fit for each AGN in Figure 14. 
The Spearman test reveals that there is no 
significant correlation between these two 
quantities, giving us confidence that our 
composite does not suffer from a bias intro- 
duced by inaccurate background subtrac- 
tion in the spectra of faint AGNs. 

7.2. Correlations of Spectral Index 
with Redshift and Luminosity 

We now turn to the question of whether 
there are any trends with redshift and/or 
luminosity in the EUV continuum slopes of 
the low-redshift AGNs in the FUSE sam- 
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pie. We show the EUV spectral slopes ver- 
sus redshift and luminosity in Figures 15 
and 16. There is a significant trend of de- 
creasing a with increasing AGN redshift 
and luminosity in the FUSE sample. The 
slopes of linear least-squares fits to a ver- 
sus log z and a versus log(AL 110 o) are listed 
in Table 4. The Spearman rank-order cor- 
relation coefficients, r s , also listed in Ta- 
ble 4, indicate that the correlations are sig- 
nificant at > 99% confidence. 

To test if these correlations are driven 
by outlier AGNs with the hardest or soft- 
est spectral slopes, we recalculated the cor- 
relation coefficients excluding the 5 AGNs 
with the softest continua from the sam- 
ple. We repeated this twice more, first 
excluding the 10 reddest AGNs, and then 
excluding the 15 reddest AGNs. We per- 
formed the same experiment with the 5, 10, 
and 15 sample AGNs with the hardest con- 
tinua. We found that the correlations of 
spectral slope with redshift and luminosity 
were robust in all cases. Examining com- 
posite spectra constructed from these same 
six subsamples, we found that although ex- 
cluding the reddest AGNs does result in a 
slightly harder composite continuum, and 
vice versa, the changes in the spectral slope 
are not significant. 

7.3. Comparisons with IUE and 
HST Samples 

7.3.1. Redshift and Luminosity Trends in 
a Combined Sample 

T02 found no significant correlation of 
O-euv with redshift or luminosity in the 
HST AGN sample. In view of the consid- 
erable correlation found in the FUSE sam- 
ple, we combine this with the sample of 



EUV spectral slopes fitted by T02 to as- 
certain whether the correlations hold for 
the combined sample. As discussed above, 
Figure 12 shows the distribution of q;euv 
in the FUSE sample and in the combined 
FUSE and HST samples. The mean and 
median spectral slopes for the FUSE and 
FUSE + HST samples are listed in Table 4. 
We use a bootstrap resampling method to 
determine the errors on the mean and me- 
dian spectral slopes, as described by T02. 
We find cwiian = -0.86 ± 0.13 for the 
FUSE sample and a mc dian = —1-40 ± 0.09 
for the combined FUSE + HST sample. 
Both values are distinctly harder than the 
medians found by T02 for the radio quiet 
and radio loud QSOs in the HST sample, 
-1.59±0.06 and -1.95±0.12, respectively. 

In Figures 17 and 18, we show the lin- 
ear fits to a versus log z and a versus 
log(ALnoo) for the FUSE sample and for 
the combined FUSE + HST sample, and 
the least-squares slopes are listed in Ta- 
ble 4. The correlations, in the sense that 
low-redshift, low-luminosity AGNs show 
harder spectral slopes, are significant at 
> 99% confidence according to the Spear- 
man rank-order test. 

7. 3. 2. Near- UV/Far- UV Spectral Break 

T02 confirmed a break at ~1200 A re- 
ported by Z97 in the spectral energy dis- 
tributions of low-redshift quasars. Their 
Figure 14 plots the individual values of 
the near-UV spectral slopes redward of 
~1200 A (ctNuv) and the far-UV spec- 
tral slopes between 500 A and the break 
(ctEuv)- This figure demonstrates that 
ajyuv is systematically larger than oceuv- 
This is also illustrated by the histograms 
of aNuv and asuv from T02 in Figure 19. 
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We seek to determine whether such a 
break is present in the low-redshift, low- 
luminosity FUSE AGN sample. The FUSE 
data do not cover the near-UV, so we 
use measurements of a^uy for our sample 
AGNs from other datasets. Many of the 
FUSE AGNs were observed with the Inter- 
national Ultraviolet Explorer (IUE) satel- 
lite. The spectral coverage of IUE, 1200- 
2000 A, is such that the IUE spectral 
slopes, ajuE, correspond roughly to a^uv 
defined by T02. Power law continuum fits 
to these spectra are compiled by S. Penton 
in an online database" 4 . 

We add to this sample of NUV spectral 
indices the NUV spectral slopes measured 
by T02 for two objects in common with 
our FUSE AGN sample, PG 1216+069 and 
PG 1543+489. To this, we also add NUV 
spectral slopes measured from combined 
STIS G140L and G230L spectra from a 
presently unpublished snapshot campaign 
to obtain HST spectra of AGNs observed 
with FUSE. There are 11 of these STIS 
spectra in common with the FUSE com- 
posite AGN sample, 8 of which were also 
observed with IUE. For the AGNs observed 
with both STIS and IUE, we use the STIS 
spectral index for our comparison with the 
FUSE data. 

In Figure 20(a), we show histograms 
of the spectral slopes in the FUSE and 
IUE/HSTb&nds for the AGNs in our com- 
parison sample. T02 found a^uv — —0.69 
versus oleuv — —1-76 for the full com- 
posite HST spectrum. The median NUV 
and EUV spectral slopes of this AGN sam- 
ple from IUE/ HST and from FUSE are 
-0.83 ± 0.04 and -0.84 ± 0.12, respec- 

4 See http: / /origins . Colorado .edu / iueagn 



tively. In Figure 20(b), we plot the spectral 
slopes we measured from the FUSE spectra 
(oiFusE) versus ansr or amE, which cor- 
respond to a nuv, to illustrate cveuv ver- 
sus cvnuv m this sample, with a line de- 
noting a fuse = c*nuv- The distribution 
of the spectral slopes about this line and 
the similar median EUV and NUV spectral 
slopes indicate that the UV spectral break 
observed in the HST composite is less pro- 
nounced or even absent in the low-redshift, 
low-luminosity FUSE AGN sample. 

While the FUSE AGN do not appear 
to show a spectral break as an ensem- 
ble, some individual objects in the sam- 
ple do show a break in the far-UV. In Fig- 
ure 21, we show two examples of AGNs 
with both FUSE and STIS G130L+G240L 
spectra, one with no break between the 
STIS and FUSE bands, and one with a 
distinct break. These two examples are 
marked as bold crosses in Figure 20(b). 
In the top panel, we show the FUSE and 
STIS spectra of PG 1322+659, a Seyfert 
1 galaxy with z = 0.168. The UV spec- 
trum of this AGN shows no evidence for a 
break, aEuv — —0.58 ± 0.29 and a^uv — 
-0.78 ± 0.04. Some low redshift AGN do 
show a spectral break, however. Kriss et al. 
(1999) reported a break in the spectrum of 
3C 273 at ~900 A from Hopkins Ultravi- 
olet Telescope data, and they interpreted 
that break as the signature of a Comp- 
tonized accretion disk spectrum (Shields 
1978; Malkan & Sargent 1982). Our com- 
parison of the FUSE and STIS spectra of 
this object, shown in Figure 21(b), confirm 
a significant break between the NUV and 
EUV spectral slopes: cxeuv — —1.60 + 0.03 
and amrv = -0.856 ± 0.005. Note that 
we simply confirm a spectral break and 
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have not attempted a self-consistent dou- 
ble power law fit to the full UV contin- 
uum here. A thorough study of the EUV- 
to-optical spectral energy distributions of 
AGNs and comparisons with accretion disk 
models using the STIS snapshot data re- 
ported above is forthcoming (Shang et al., 
in preparation). 

7.4. Correlation of Spectral Index 
with Accretion Properties 

Here we investigate whether the spec- 
tral slopes of the low-redshift AGNs in the 
FUSE sample are correlated with observ- 
able properties of the central black hole or 
the mass accretion process. For 22 AGNs 
in the FUSE sample, we have estimates of 
black hole mass from Kaspi et al. (2000) 
and M c Lure & Dunlop (2001). This quan- 
tity is listed in column 6 of Table 1 and is 
plotted against the spectral indices of indi- 
vidual AGNs in Figure 22. We find a cor- 
relation of spectral slope with black hole 
mass, at 96% confidence according to the 
Spearman rank-order test. The trend runs 
in the sense that AGNs with lower black 
hole masses show harder spectral slopes. 
A linear least-squares fit gives a slope of 
—0.67 ± 0.35, shown by the line in Fig- 
ure 22. 

For standard, geometrically thin, opti- 
cally thick accretion disks, the tempera- 
ture of the disk is proportional to its lu- 
minosity and to the mass of the central 
black hole as: T ~ L 1 / 4 M B ^' 2 (Shakura 
& Sunyaev 1973). Therefore, the thermal 
disk emission will peak further in the blue 
in the spectra of AGNs with lower mass 
black holes. Thus, the spectral index ob- 
served in the NUV will extend through the 
EUV, and we may expect a correlation be- 



tween a and temperature. We show a ver- 
sus L^M^J 2 in Fi gure 23 and the linear 
least-squares fit with a slope 0.86 ± 0.91. 
The slope has the expected sign, but the 
correlation is not statistically significant, 
70%. 

8. Discussion 

The FUSE AGN sample is approxi- 
mately one-half the size of the HST sample 
presented by T02, and it is distinctly differ- 
ent in terms of both redshift (-2 mc dian = 0.10 
versus z median = 0.96 for the HST AGNs) 
and luminosity (median log(AL 110 o) = 45.0 
versus median log(ALnoo) = 45.9). Like 
the Z97 and T02 HST samples, this FUSE 
archival sample is a heterogeneous sam- 
ple of AGNs observed with FUSE for var- 
ious reasons. Presumably, however, most 
AGNs in the sample were observed because 
they were known a priori to be bright in 
the UV, particularly in the NUV. How- 
ever, we note that Figures 5 and 6 illus- 
trate that the sample is well-populated 
with intrinsically low-luminosity AGNs. 
The EUV spectral index of the FUSE 
composite is significantly harder than the 
HST composite, a = — 0.56^o;<l versus 
a = -1.76 ±0.12. 

In constructing our composite spec- 
trum, we followed the methodology of Z97 
and T02 by using uniform weighting for 
each AGN spectrum. This method pre- 
vents the overall composite spectrum from 
being dominated by the brightest AGN 
with spectra of the highest S/N, such as 
3C273 (a = —1.60). This in turn provides 
the best estimate of the spectrum of the 
average AGN (Z97) and the best compari- 
son with the HST results of Z97 and T02. 

At face value, the harder spectral slope 
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of a composite formed from low-redshift 
AGNs seems to lend credence to a scenario 
in which the spectral break in the HST 
composite spectrum reported by Z97 and 
T02 is caused by an undercorrection for in- 
tergalactic absorption, as proposed by Bi- 
nette et al. (2003). However, this scenario 
does not easily explain the correlation be- 
tween spectral slope and AGN luminosity 
we report, nor does it account for the en- 
hanced emission lines in the FUSE com- 
posite relative to the HST composite. 

A more natural explanation of these re- 
sults is that we are seeing a manifesta- 
tion of the physical mechanism put forth 
to explain the Baldwin effect. This corre- 
lation is generally attributed to the ten- 
dency for low-luminosity AGNs tend to 
show harder ionizing continua (Zheng & 
Malkan 1993; Wang et al. 1998; Diet- 
rich et al. 2002). This suggests that both 
the enhanced high-ionization emission line 
strengths and the harder continuum shape 
of the FUSE composite spectrum with re- 
spect to the HST composite are due to the 
larger fraction of relatively low-luminosity 
AGNs in the FUSE sample. T02 explain 
the excess of C IV emission in the HST 
composite relative to the SDSS composite 
(Vanden Berk et al. 2001) in the same way. 
The larger equivalent widths of the high- 
ionization emission lines O VI and Ne Vlll 
in the FUSE composite relative to the HST 
composite is the best evidence that the 
harder spectral index of the FUSE com- 
posite is not caused by any wavelength- 
dependent systematic error in our correc- 
tions for Galactic and intergalactic absorp- 
tion. In addition, the significant corre- 
lation between individual AGN spectral 
slopes and luminosities for both the FUSE 



sample and the combined FUSE + HST 
samples supports the Baldwin effect inter- 
pretation. 

Z97 interpreted the spectral break at 
1050 A in their HST composite, along with 
a depression of the flux near the Lyman 
limit, as the signature of a Comptonized 
accretion disk spectrum. Relativistic ef- 
fects on the emergent spectrum from an 
inclined accretion disk as well as Comp- 
tonization from a hot corona can broaden 
an intrinsic Lyman limit break (Lee, Kriss, 
& Davidsen 1992). Z97 estimated a Ly- 
man limit optical depth r = 0.8 for a 10% 
depression in the flux. We have shown 
that the AGNs in this low-redshift, low- 
luminosity FUSE sample tend to lack such 
a spectral break, and we find no evidence 
for a depression in the spectrum blueward 
of the Lyman limit in the FUSE com- 
posite. From the S/N in the compos- 
ite, we estimate that any spectral discon- 
tinuity is less than 10%. Using non-LTE 
accretion disk models with no Compton 
scattering, Hubeny et al. (2000) demon- 
strated that emission may wash out Ly- 
man edge discontinuities in QSOs with 
L > O.OlLEdd) particularly in systems with 
M BH < 1O 9 M . We calculate the Edding- 
ton ratio for the 22 FUSE AGNs for which 
we have estimates of black hole mass in 
Table 1 using that black hole mass and 
bolometric luminosities from Woo & Urry 
(2002) or Padovani & Rafanelli (1988), also 
listed in Table 1. The Eddington ratio 
is plotted versus black hole mass in Fig- 
ure 24, which indicates that the FUSE 
AGN occupy the region of parameter space 
noted above. The lack of a flux disconti- 
nuity in the FUSE composite is consistent 
with the Hubeny et al. (2000) models in 
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this sense. However, their models do pre- 
dict a change in the spectral slope for these 
AGNs at wavelengths blueward of the Ly- 
man limit, a break we do not find in the 
FUSE composite. 

The EUV spectral slope found here for 
AGNs with z < 0.67 does imply a spectral 
break between the far UV and the soft X- 
rays given the soft X-ray slopes reported 
by Laor et al. (1997), who find < a x >= 
-1.72 ± 0.09 and < a x >= -1.15 ± 0.14 
for radio-quiet and radio-loud low-redshift 
QSOs, respectively. It is consistent with 
the presence of a soft X-ray excess and 
steep soft X-ray power law to join with 
the 0.2-2 keV power law, as Mathews & 
Ferland (1987) inferred from the observed 
strength of He II A 1640 emission. The hard 
slope of the FUSE composite may also help 
resolve questions raised by Korista, Fer- 
land, & Baldwin (1997) about how a soft 
EUV continuum such as the one found by 
Z97 and T02 could account for the ob- 
served strength of this line. Neglecting any 
turnover in the EUV-to-soft-X-ray spectral 
energy distribution, the FUSE composite 
predicts about eight times as many pho- 
tons at 4 Ryd as the a = —2 power law 
these authors used to calculate He II A1640 
equivalent widths of 0.6-0.8 A for cover- 
ing fractions of 10%. The FUSE data do 
not cover rest frame 1640 A, and although 
He II A1085 emission line may be present in 
the emission line complex identified in the 
composite at 1084 A (see § 6), this iden- 
tification is uncertain. We therefore leave 
a more complete consideration of the im- 
plications of the FUSE EUV spectral slope 
on AGN emission lines for future work. 

If AGNs build up black hole mass via 
accretion over their lifetimes, the evolu- 



tionary models of Wandel (1999a,b) pre- 
dict higher mass (luminosity) AGNs will 
have lower temperature accretion disks and 
softer spectral slopes in the UV and X- 
rays. The absence of a spectral break in 
the FUSE composite can be explained by 
such models. Low-luminosity AGNs pos- 
sess hotter accretion disks, and this shifts 
the UV bump and the Compton break to 
shorter wavelengths (Z97). The anticorre- 
lation we find between EUV spectral slope 
and black hole mass estimates compiled 
from the literature also supports this in- 
terpretation. This is undermined some- 
what by the low significance of the cor- 
relation between EUV spectral index and 
accretion disk temperature, estimated by 
L 1 / 4 M B 1 J 2 . For AGNs accreting at the Ed- 
dington limit, L oc M BH and T oc M^ 1 / 4 
and the significant correlation between a 
and M BH may directly reflect a trend of a 
and accretion disk temperature. The Ed- 
dington ratios listed in Table 1 and plotted 
in Figure 24 indicate that in fact most of 
the AGNs in the sample shine at luminosi- 
ties below the Eddington limit. However, 
we note that the estimate we use for the 
accretion disk temperature is a combina- 
tion of two observables. If a trend between 
spectral slope and disk temperature is in- 
deed present, it could easily be swamped 
by the considerable uncertainties in the lu- 
minosity and black hole mass. 

Finally, we discuss the dip in the com- 
posite spectrum blueward of the Ne vm 
emission line. Possibly, this arises from the 
superposition of absorption features aris- 
ing from highly ionized gas along the line of 
sight, either intrinsic to the AGNs or in the 
IGM. It is unlikely to reside in intergalactic 
absorbers comprising the warm-hot IGM 
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(WHIM). These absorbers have been ob- 
served in VI in several QSO sightlines 
(Tripp, Savage, & Jenkins 2000; Oegerle 
et al. 2000; Tripp & Savage 2000; Tripp 
et al. 2001; Sembach et al. 2001; Sav- 
age et al. 2002). We see no correspond- 
ing dip in the flux at the position of the 
O VI doublet in the FUSE composite spec- 
trum. Whether these absorbers are pho- 
toionized or collisionally ionized, it is likely 
that O vi/Ne vill > 1 (Tripp & Savage 
2000). On the other hand, if the absorb- 
ing gas is collisionally ionized and 5.75 < 
logT < 6.00, or if nonequilibrium condi- 
tions hold, the abundance of Ne Vill could 
be appreciable, even comparable to that of 
O VI (Heckman et al. 2002). At such tem- 
peratures, O VII and even O Vill absorp- 
tion would be expected and have been ob- 
served (Nicastro et al. 2002; Mathur, Wein- 
berg, & Chen 2003). Nonetheless, if this 
flux depression is attributable to Ne Vill 
in the WHIM, the question of why we find 
no depression due to intervening O VI is 
still open. If we do attribute the dip to 
intervening Ne Vill absorption, we may es- 
timate the mean cosmological mass density 
in these absorbers using 



tt b (Ne vm) = 
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(Tripp et al. 2000), where /i is the mean 
atomic weight, taken to be 1.3, /(Ne vm) 
is the ionization fraction of Ne vm, and 
(Ne/H) is the neon abundance by number. 
The quantity ^jiVj(Ne vm) is the total 
column density of Ne Vill estimated from 
the apparent column density in the absorp- 
tion dip in the composite. The effective 
optical depth of the feature implies a total 



column density N(Ne vm) ~ 6 x 10 15 cm 2 , 
assuming optically thin absorption (Sav- 
age & Sembach 1991). The term AX is 
the total absorption distance probed by all 
the sample AGN assuming q = 0.5 (Bah- 
call & Peebles 1969). For [Ne/H]=-1 and 
/(Ne Vlll)=0.2, the peak fractional abun- 
dance of Ne vm at logT = 5.85 (Shapiro & 
Moore 1976), we estimate ~ 0.027/1^5. 
This value is a factor of ~10 larger than 
the more careful lower limits set from O VI 
absorbers (Tripp & Savage 2000; Savage et 
al. 2002), and it would imply that over 60% 
of the baryons at z < 0.7 reside in Ne vm 
absorbers. 

This flux depression has a symmet- 
ric appearance with a velocity extent of 
~10,000 km s -1 , and the centroid of the 
absorption lies within ~17,000 km s _1 of 
the Ne Vill emission line. It is difficult 
to imagine a redshift distribution of in- 
tervening absorbers that would naturally 
give rise to these characteristics, which are 
more closely akin to broad-absorption-line 
(BAL) systems. It is more likely that this 
feature is associated with the AGNs them- 
selves. Broad, intrinsic highly- ionized ab- 
sorption from Ne Vill has been observed 
in several high redshift QSOs (Q 0226- 
1024: Korista et al. 1992; SBS 1542+541: 
Telfer et al. 1998; and PG 0946+301: Arav 
et al. 1999), as have narrow associated 
absorption features (UM 675: Hamann 
et al. 1995, 1997; HS 1700+6416: Petit- 
jean, Riediger, & Rauch 1996; J2233-606: 
Petitjean & Srianand 1999; and 3C 288: 
Hamann, Netzer & Shields 2000). We 
find no obvious Ne VIII absorption fea- 
tures in any single AGN spectrum. We 
deliberately excluded AGNs known to be 
BALs from the sample, but many AGN 
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showing narrow intrinsic UV absorption 
troughs associated with X-ray "warm ab- 
sorbers", such as Mrk 279 (Scott et al. 
2004), Mrk 509 (Kriss et al. 2000a; Krae- 
mer et al. 2003), NGC 3783 (Gabel et 
al. 2003), and NGC 7469 (Kriss et al. 
2000b; Kriss et al. 2003) are included in 
the FUSE composite sample. In photoion- 
ization models of intrinsic absorbers, Ne 7+ 
and 6+ can coexist over large regions 
(Hamann 1997; Wang et al. 2000). Fur- 
thermore, given solar abundances and the 
similar oscillator strengths of these two 
doublets, we might expect to find a corre- 
sponding dip in the composite from associ- 
ated O VI absorption. However, the BAL 
in the spectrum of SBS 1542+541 shows 
absorption due to highly-ionized species, 
including Ne vill, and lacks absorption fea- 
tures from low-ionization species such as 
C ill and Si IV seen in other BALs (Telfer 
et al. 1998). Their results suggest that 
the covering fractions of various ions can 
be dependent upon creation ionization po- 
tential. They do find O VI absorption in 
this BAL system, but our selection against 
known BALs in the FUSE sample from fea- 
tures like C IV and O VI, combined with 
an ionization-dependent covering fraction 
is a plausible explanation of the flux de- 
pression blue ward of Ne VI 1 1 in the FUSE 
composite, and the lack of one blueward of 
O vi. 

9. Summary 

We summarize our results as follows: 

1. We construct a composite EUV (630- 
1155 A) spectrum of AGNs with z < 
0.67 from archival FUSE data. 



2. We find that the best-fit spectral in- 
dex of the composite is a = — 0.56^Q'j§ 
The conservative estimate of the to- 
tal error in the spectral index in- 
cludes the standard deviation in a in 
1000 bootstrap samples of the FUSE 
data set, uncertainties in the extinc- 
tion correction applied to the FUSE 
spectra, and uncertainties in the col- 
umn density distribution parameter 
of the intervening Lya forest. 

3. We find that O Vl/Ly/3 and Ne vm 
emission are enhanced in the FUSE 
composite relative to the HST com- 
posite. 

4. The FUSE composite is significantly 
harder than the EUV portion of the 
HST composite spectrum of T02 who 
find a = —1.76 ±0.12 for 332 spectra 
of 184 AGNs with z > 0.33. 

5. We find significant correlations of 
EUV spectral index with redshift and 
luminosity in the FUSE AGN sample 
and in the FUSE sample combined 
with the sample of T02 in the sense 
that lower redshift/luminosity AGNs 
have harder spectral slopes. 

6. From comparisons of the FUSE spec- 
tra with data from HST and IUE, we 
find no evidence for a UV spectral 
break at ~1200 A in the FUSE AGN 
sample. 

7. We find a significant correlation of 
EUV spectral index with black hole 
mass for 22 AGNs for which such 
measurements exist in the litera- 
ture. The trend of EUV spectral 
index with disk temperature, esti- 
mated from the AGN luminosities 
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and black hole masses, runs in the 
expected direction. The correlation 
between these quantities is not signif- 
icant, perhaps due to observational 
uncertainties. 

8. The FUSE sample is dominated by 
low-redshift, low-luminosity AGNs. 
Therefore, items 3-7 suggest the 
physical underpinnings of the Bald- 
win effect, the tendency for lower 
luminosity AGNs to have hotter ac- 
cretion disks and harder spectra in 
the UV and soft X-rays. 
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Table 1 
AGN in FUSE Sample 



Name 


z 


E(B-V) 


log(AL) 1 


«a 2 


F A 2 


logM BH /M 3 


log(L bol ) 4 


NGC 3783 


0.010 


0.119 


43.7 


1.30 ± 0.14 


2.04 ± 0.03 


7-041^9 


44.41 


Mrk 352 


0.015 


0.061 


43.1 


2.69 ± 0.42 


2.37 ± 0.01 






NGC 7469 


0.016 


0.069 


43.9 


1.76 ± 0.09 


1.32 ± 0.01 


r oo+O 43 

6.881^3 


45.28 


Mrk 79 


0.022 


0.071 


43.6 


1.17 ± 0.41 


0.29 ± 0.01 


o m +0 17 
».U1_ 24 


44.57 


Ark 564 


0.025 


0.060 


43.3 


1.57 ± 0.18 


0.160 ± 0.004 






Mrk 335 


0.026 


0.035 


44.3 


1.68 ± 0.06 


1.20 ± 0.01 


n ro4-0 16 


44.69 


Mrk 290 


0.030 


0.015 


43.5 


1.90 ± 0.23 


0.14 ± 0.04 






Mrk 279 


0.030 


0.016 


44.5 


1.91 ± 0.05 


1.270 ± 0.007 






Mrk 817 


0.031 


0.007 


44.4 


1.50 ± 0.08 


0.87 ± 0.01 


n r- 1 3 


44.99 


Ark 120 


0.032 


0.128 


44.5 


2.82 ± 0.14 


1.25 ± 0.02 


°- z '-o.w 


44.91 


IRAS F11431-1810 


0.033 


0.039 


44.4 


4.03 ± 0.18 


1.13 ± 0.02 






Mrk 509 


0.034 


0.057 


44.4 


1.04 ± 0.28 


1.34 ± 0.05 




45.03 


Mrk 618 


0.035 


0.076 


44.5 


2.51 ± 0.34 


1.20 ± 0.06 






ESO 141-G55 


0.036 


0.111 


44.8 


3.27 ± 0.14 


2.78 ± 0.06 


8.851^3 


45. 62 4 


Mrk 9 


0.040 


0.059 


44.3 


0.95 ± 0.18 


0.344 ± 0.009 






1H 0707-495 


0.041 


0.095 


44.4 


3.63 ± 0.11 


0.86 ± 0.01 






NGC 985 


0.042 


0.033 


44.2 


3.10 ± 0.06 


0.403 ± 0.003 






KUG 1031+398 


0.042 


0.015 


43.3 


2.83 ±0.27 


0.051 ±0.001 






Mrk 506 


0.043 


0.031 


44.0 


2.42 ±0.24 


0.209 ±0.007 






Fairall 9 


0.047 


0.027 


43.9 


1.04 ±0.33 


0.172 ±0.008 


7 Q9+ - 13 

' - yz -0.22 


45.23 


Mrk 734 


0.050 


0.032 


44.5 


1.37±0.19 


0.54 ±0.01 






IRAS 09149-62 


0.057 


0.182 


45.0 


4.34 ±0.32 


1.31 ±0.05 






ESO 265-G23 


0.056 


0.096 


44.4 


2.14 ±0.43 


0.38 ±0.01 






3C 382 


0.058 


0.070 


44.4 


2.48 ±0.29 


0.36 ±0.01 






PG 1011-040 


0.058 


0.037 


44.5 


0.67 ±0.07 


0.401 ±0.003 







Table 1 — Continued 



1\Tp TY1P 


As 


EfR-Vl 






Ex 2 


lug, ivibh/ iviq lug^boU 


Mrk 1298 


0.060 


0.055 


44.5 


1.14 + 0.46 


0.39 + 0.02 




I Zw 1 


0.061 


0.065 


44.2 


-0.51 + 0.27 


0.157 + 0.006 




Ton S180 


0.062 


0.014 


44.9 


1.93 + 0.14 


0.96 + 0.01 




II Zw 136 


0.063 


0.044 


44.6 


1.57 + 0.31 


0.464 + 0.006 




PG 1229+204 


0.063 


0.027 


44.4 


-1.12 + 0.58 


0.278 + 0.007 


7.9318$ 45 - 01 


Ton 951 


0.064 


0.037 


44.8 


0.36 + 0.17 


0.720 + 0.005 




MR 2251-178 


0.066 


0.039 


44.5 


0.69 + 0.08 


0.299 + 0.003 




Ton 1187 


0.070 


0.011 


44.4 


0.84 + 0.59 


0.223 + 0.006 




Mrk 205 


0.071 


0.042 


44.2 


1.76 + 0.60 


0.107 + 0.003 




Mrk 478 


0.079 


0.014 


44.8 


1.29 + 0.32 


0.443 + 0.005 




VII Zw 118 


0.080 


0.038 


44.8 


1.46 + 0.30 


0.427 + 0.004 




PG 1211+143 


0.081 


0.035 


45.2 


2.83 + 0.27 


0.868 + 0.008 


7.371^3 45.81 


Mrk 1383 


0.086 


0.032 


45.2 


0.72 + 0.08 


1.06 + 0.004 




PG 0804+761 


0.100 


0.035 


45.7 


-0.53 + 0.05 


1.010 + 0.003 


8.2118$ 45.93 


PG 1415+451 


0.114 


0.009 


44.7 


1.88 + 0.24 


0.155 + 0.003 




Ton S210 


0.116 


0.017 


45.4 


-0.57 + 0.07 


0.790 + 0.007 




RX J1230.8+0115 


0.117 


0.019 


45.3 


-0.24 + 0.20 


0.63 + 0.01 




Mrk 106 


0.123 


0.028 


45.0 


-0.16 + 0.22 


0.280 + 0.005 




Mrk 876 


0.129 


0.027 


45.7 


-0.92 + 0.05 


1.070 + 0.005 




PG 1626+554 


0.133 


0.006 


45.0 


-0.04 + 0.23 


0.218 + 0.004 




Q0045+3926 


0.134 


0.052 


45.0 


1.10 + 0.31 


0.210 + 0.005 




PKS 0558-504 


0.137 


0.044 


45.5 


1.14 + 0.06 


0.695 + 0.003 




PG 0026+129 


0.142 


0.071 


45.3 


1.24 + 0.16 


0.445 + 0.006 


7.42l°;°9 45.39 


PG 1114+445 


0.144 


0.016 


44.3 


3.44 + 0.67 


0.033 + 0.002 




PG 1352+183 


0.152 


0.019 


44.7 


1.89 + 0.40 


0.080 + 0.002 





Table 1 — Continued 



Name 




z 


E(B-V) 


log(AL) 1 




«A 2 








Fa 


2 




logM BH /M 3 


log(I 


boi) 4 


MS 07007+6338 





153 





051 


45 


3 





68 


± 





16 


0.429 


± 





004 








PG 1115+407 





154 





016 


44 


8 





89 


± 





23 


0.099 


± 





002 








PG 0052+251 





155 





047 


45 


2 





90 


± 





03 


0.315 


l 





005 


,,0+0.13 

8.481^ 


/i r 

45 


93 


T>/"^ i on^ i nor 

PG 1307+085 





155 





034 


45 


5 


/A 




07 


± 





08 


0.557 


1 

± 


/A 




r\r\ a 

004 


o co+0.16 
°- dz -0.29 


a r 

45 


83 


30 273 





158 





021 


46 


4 





40 


± 





03 


3.630 


l 

± 


/A 




007 


o cr /I +0.41 


/I T 

47 


o A 4 

24 


PG 1402+261 





164 





016 


45 


5 





82 


± 





13 


0.431 


1 

± 





r\r\ a 

004 


« ~ r +0.91 3 

7- 76^ Jo 


a r 

45 


13 


PG 1048+342 





167 





023 


45 





1 


18 


± 





26 


0.138 


1 





003 








T /' I -1 f ) ,A ,A . Z" 1 P" /A, 

PG 1322+659 





168 





019 


45 


1 


1 


42 


± 





29 


0.177 


l 


(A 




r\r\ a 

004 








PG 2349-014 





-1 T ,1 

174 





027 


45 


4 


— 1 


08 


± 





1 A 

14 


0.258 


l 


(A 




/A/A O 

003 


n n/^+0.91 3 

9-26-0.30 


/I V 

45 


94 


PG 1H6+215 





176 





023 


45 


9 


— 


40 


± 





03 


0.891 


1 

± 


(A 




003 


f 7ri _|_o.91 3 
8 - 70 -0.30 


46 


/• > , A 

02 


T ' T ~) S \ C~i Tf\-i r r /A /A ,A ,A 

FBQS J2155-0922 





190 





046 


46 








02 


± 





12 


1.03 


± 





01 








4C +34.47 





206 





037 


45 


1 





12 


± 





36 


0.099 


± 





003 








ta s~\ r\r\ a r~7 , t) (\/ ■ 

PG 0947+396 





OAP 

206 





019 


45 


2 


1 


15 


± 


/A 




25 


0.155 


1 


/A 




AAO 

003 








TTT — \ 1 ATA n^71 1 

HE 1050-2711 


/A 




208 





070 


45 


2 


/A 




55 


± 





30 


0.146 


l 





004 








TTTH 1 1 1 r 1 *"70 r - 

HE 1115-1735 





217 





037 


45 


6 


— 


16 


± 





"1 o 

12 


0.278 


1 





003 








PG 0953+414 





O O /I 

234 





013 


46 





/A 




75 


± 





03 


/A pnr 

0.605 


1 





001 


o qi +0.07 
°- Zi -0.08 


a n 

46 


16 


TTT7 1 1 ni r i pi o 

HE 1015-1618 





O /I T 

247 





078 


45 


4 


— 1 


33 


± 





34 


0.135 


1 





r\r\ a 

004 








t ^ / i 1/1/1/1 i ,1 /at 

PG 1444+407 





267 





014 


45 


6 





91 


± 


/A 




13 


0.202 


1 

± 





002 


o cr a +0.91 3 

8.54I 42 


A r 

45 


93 


PKS 1302-102 





278 





043 


45 


8 





09 


± 





05 


0.280 


± 





001 


- n +f| Q4 o 

°- 'O-0.30 


45 


86 


HE 0450-2958 





286 





015 


45 


5 





79 


± 





13 


0.157 


± 





002 








PG 1 100+772 





311 





034 


45 


7 





22 


± 





15 


0.238 


± 





005 








Ton 28 





330 





022 


45 


9 


1 


97 


± 





13 


0.307 


± 





005 








FBQS J083535.8+24 





331 





031 


45 


7 


2 


42 


± 





25 


0.188 


± 





005 








PG 1216+069 





331 





022 


45 


7 


1 


20 


± 





18 


0.184 


± 





004 








PG 1512+370 





371 





022 


45 


6 


2 


77 


± 





53 


0.137 


± 





009 









Table 1 — Continued 



l\Ta TTIP 
1 >l dill" 




y 

A* 


1 J \ I J V J 




m 2 




Fx 2 


lug 1V1BH/ lvi iUg^bol 


PG 1543+489 





400 


0.018 


45.7 


0.83 + 


54 


0.099 + 


007 


FBQS J152840.6+28 





450 


0.024 


46.0 


0.61 + 


19 


0.177 + 


004 


PG 1259+593 





478 


0.008 


46.1 


0.42 ± 


02 


0.1810 + 


0004 


HE 0226-4110 





495 


0.016 


46.4 


0.70 ± 


04 


0.345 ± 


001 


HS 1102+3441 





510 


0.024 


46.0 


2.83 + 


24 


0.154 + 


003 


RX J2241.8-44 





545 


0.011 


46.0 


3.47 + 


16 


0.122 + 


002 


PKS 0405-12 





573 


0.058 


46.9 


3.35 + 


08 


1.13 + 


01 


HE 1326-0516 





578 


0.030 


46.4 


-0.33 + 


13 


0.141 + 


002 


HE 0238-1904 





631 


0.032 


46.5 


1.61 + 


11 


0.287 + 


003 


3C 57 





669 


0.022 


46.4 


1.52 + 


14 


0.169 + 


003 



1 uminosity at 1100 A, log(ALnoo) for fi = 1 and H = 60 km s 1 Mpc 1 , see Equ. 2 

2 Spectral index, F x oc \~ ax , and flux at 1000 A in units of 10~ 13 ergs s" 1 cm -2 A -1 from power law 
continuum fit 

Measurements of M B h from Kaspi et al. (2000) unless marked with "3", which are from M c Lure & Dunlop 
(2001) 

4 Bolometric Luminosity calculated by Woo & Urry 2002 unless marked with "4" , which are from Padovani 
& Rafanelli (1988) 



Table 2 

Fits to Full Sample and Subsample 
Composites 



Sample 


No. Spectra 


No. AGNs 


a 


Full sample 


128 


85 


-0.56 


z < 0.10 


65 


39 


-0.74 


z > 0.10 


63 


46 


-0.74 


log(ALnoo) < 45.0 


64 


40 


0.19 


log(AL n oo) > 45.0 


64 


45 


-0.84 


E(B - V) < 0.032 


65 


44 


-0.50 


E(B-V) > 0.032 


63 


41 


-0.21 



Table 3 
Emission Lines 



Line 


Ao (A) 


Flux 


EW (A) 


Nevm + OlV 


772 


6.3 ±0.3 


8.4 


Om 


831 


1.6 ±0.3 


2.4 


H I Ly series + Svi 


944 


3.2 ±0.5 


5.6 


Ly7 


973 


1.1 ±1.1 


2.0 


Cm 


976 


1.5 ±0.4 


2.9 


Nm 


991 


1.1 ±0.2 


2.1 


Ly/3 


1026 


6.2 ±1.3 


12.6 


Ovi 


1033 


10.0±0.3 


20.3 


Siv 


1062 


1.7 ±0.3 


3.6 


Siv 


1073 


0.8 ±0.1 


1.8 


Nil + Hen + An 


1084 


0.23±0.05 


0.4 



23 



Table 4 
Fits to Individual AGN 











a vs. log 




a vs. log 


■\L 


Sample 


N 


Omean 


^med. 


Slope 


r s 


Slope 


r s 


FUSE 

FUSE + HST 


85 
164 


-0.74±0.13 
-1.25±0.09 


-0.86±0.13 
-1.40±0.09 


-0.83±0.02 
-0.88±0.13 


-0.31 
-0.45 


-0.56±0.01 
-0.49±0.09 


-0.36 
-0.38 
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Fig. 1. — Number of spectra used in the composite spectrum as a function of rest wavelength. 
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Fig. 2. — Signal-to-noise ratio in the composite spectrum per 1 A bin. 
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Fig. 3. — Top panel: Composite AGN spectrum with power law continuum fit {dashed lines) 
and wavelength regions used in fit {solid line segments); HST composite from Telfer et al. 
(2002) shown for comparison. Bottom panel: Ratio of FUSE to HST composite spectra. 
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Fig. 4. — Lyman valley correction for a z — 0.1 AGN with the parameters used by T02, and 
the correction used here, where 7 = 0.15, /3 — 2.0 in Equ. 1. The dashed lines show how the 
correction changes for different 7 and (3. 
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Fig. 5. — Luminosity versus redshift for FUSE and HST AGN with lines marking median 
redshift (z = 0.10) and median luminosity (logAL 1100 = 45.0). 
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Fig. 6. — Histograms of redshifts and luminosities of FUSE and HST AGN with vertical 
lines marking median FUSE redshift and luminosity. 
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Fig. 7. — Composites made from spectra of AGNs with z < 0.10 and z > 0.10 with power 
law continuum fits (dashed lines) and wavelength regions used in fit (solid line segments). 
Overall composite with power law continuum fit shown for comparison. 
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Fig. 8. — Composites made from spectra of AGNs with log(ALnoo) < 45.0 and log(ALnoo) > 
45.0 with power law continuum fits (dashed lines) and wavelength regions used in fit (solid 
line segments). Overall composite with power law continuum fit shown for comparison. 
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Fig. 9. — Histogram of E(B — V) values for FUSE AGNs. The solid line marks the median 
E(B-V), 0.032. 
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Fig. 10.— Distribution in redshift and 1100- A luminosity for AGNs with E(B - V) > 0.032 
(filled Mangles) and AGNs with E(B - V) < 0.032 (crosses) 
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Fig. 11. — Composite AGN spectrum with fits to continuum {dashed line) and continuum + 
emission lines (solid line). See Table 3 for fit parameters. 
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Fig. 12. — Histogram of EUV spectral slopes in FUSE sample and in the combined FUSE 
and HST sample. 
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Fig. 13. — EUV spectral index versus E(B — V) for FUSE AGNs with best linear least- 
squares fit. 
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Fig. 14.— EUV spectral index versus flux at 1000 A for FUSE AGNs. 
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Fig. 15. — EUV spectral index versus redshift with best linear least-squares fit. 
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Fig. 16. — EUV spectral index versus 1100- A luminosity with best linear least-squares fit. 
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Fig. 17. — Same as Figure 15 for both FUSE (solid squares) and HST (Radio-loud: solid 
circles, Radio-quiet: open squares) samples with best linear least-squares fit to FUSE sample 
(solid line) and best linear least-squares fit to combined FUSE + HST sample (dashed line). 
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Fig. 18. — Same as Figure 16 for both FUSE (solid squares) and HST (Radio-loud: solid 
circles, Radio-quiet: open squares) samples with best linear least-squares fit to FUSE sample 
(solid line) and best linear least-squares fit to combined FUSE + HST sample (dashed line). 
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Fig. 20. — (a) Histogram of spectral slopes measured with IUE for rest wavelengths 1200- 
2000 A or EST for rest wavelengths 1150-3200 A (olnvv) an d those measured with FUSE 
for 900 < A < 1200 A ((Xfuse)- (b) Comparison of cxfuse and a^uv- Filled squares 
denote a^i/v's measured from IUE spectra, open squares denote o/vuv's measured from 
EST spectra, and bold crosses mark the two AGNs shown in Figure 21. Dotted line shows 

O-FUSE = 0>NUV- 
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Fig. 21. — Comparison of HST spectra {light grey) and smoothed FUSE spectra {dark grey) 
of PG 1322+659 {top panel) and 3C 273 {bottom panel). The solid lines show the power law 
continuum fits, and the dotted lines indicate the la uncertainties on the power law index. 
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Fig. 22. — EUV spectral index versus black hole mass listed in Table 1 for 21 AGNs in the 
FUSE sample, with best linear least-squares fit. 
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Fig. 23. — EUV spectral index versus temperature (~ log(L 1 / 4 M B ^/ 2 )) with best linear 
least-squares fit. 
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Fig. 24. — Eddington ratio versus black hole mass for 20 FUSE AGNs. See Table 1. 
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